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The exponentially-mapped-past approach 


It is often desirable to compute something like an arithmetic mean 





for a (time) continuous variable x(t). A very simple approach could look like this: 


ty 


= 1 
r= f o(t)at 


to 





Although this approach works fine on an analog computer it requires fixed times to 
and t, which is only practical in few cases. 

To overcome this problem [OTTERMAN(1960)] introduced exponentially-mapped- 
past (EMP) variables in order to extend the idea of a mean to continuous variables in 
continuous time on which the classic EAI application note [EAI 1.3.2(1964)] is based.* 
The basic idea is to introduce a weighting function that ensures that recent values 
influence the result more strongly than values in the past. The following equation 
demonstrates this technique with the integral running from the most distant past —oo 
to 0 (=“now’):? 


0 
z(0) =a / a(t)e™ dt 
Here, a denotes a normalization factor: 
; 1 
/ eS — 
a 





1OTTERMAN’s work has its roots in [FANO(1950)]. This application note closely follows these sources as well as 
[EAI 1.3.2(1964)]. 
2% denotes the EMP mean. 
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Figure 1: EMP mean circuit — the parameter a determines how quickly the weigthing function 
“forgets” past input values 


Implementing this scheme on an analog computer is straightforward and guarantees 
that the integrator will not overload even with long run times (given that x(t) remains 
in suitable bounds). 

This can be further generalized as 


T T 
z(T) =a ; a(t)e*’—? dt = ae°T / x(t)e™ dt, (1) 


a convolution integral of the input function x(t) and an exponentially decaying weigth- 
ing function, the derivative of which with respect to 7’ is 


iP 
ai =a | —ae°? / a(t)e dt +e eT (T) | = ax(T)—az(T). (2) 


Based on (2) the analog computer setup shown in figure 1 can be directly derived. 
This basically implements a “leaky integrator” which can also be seen as a low-pass 
RC filter. It should be noted that this only works if no exact estimation of the mean 
value is required during the startup time of the computation. After a step input, the 
output will reach 95% of the step height in the time interval 3/a. This must be taken 
into account for the startup time as well. 
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Figure 2: EMP variance circuit 

















One can now ask if the idea of EMP variables also allows the calculation of a type 
of variance. In the discrete case the variance is defined by 





Analogously to (1) this can be extended to continuous variables in continuous time as 


o2(T) =a / (a(t) — &(t))? eX dt = ae“°? / (x(t) — &(t))? e dt 


which can be mechanized by the analog computer setup shown in figure 2. 
Equally straightforward is the computation of an EMP autocorrelation p(7) based 


on 
T 


pt) =e / x(t)a(t — pew 4 dt 
as shown in figure 3 where 7 represents the time delay used for the correlation. The 
time delay function shown can be implemented using various techniques such as PADE 
or STUBBS-SINGLE approximations (cf. [ULMANN(2020), pp. 95-108]). 
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Figure 3: EMP autocorrelation circuit 


The WIENER-KHINCHIN theorem states that the spectral decomposition of the 
autocorrelation function of a suitable function is given by the power spectrum of that 
function.? Thus, it is possible to compute an EMP power spectrum based on (rT). 

The EMP FOURIER transform is defined as 


T TE 
F(w) —a / rte OF Dei dt = qe it / a(t)e oT) giv (F—-4) dt. 


The power spectrum is P(w) = |F'(w)|?, i.e. in the EMP case it is 


T 2 


fewere sin(w(T' — t)) dt 


(oe) 


The corresponding analog computer setup is shown in figure 4. At its heart is a simple 
quadrature oscillator consisting of two integrators and a summer in a loop. This yields 
both, the sine and cosine components, the squares of which are summed to yield the 
desired output. 





3Questions regarding convergence criteria are beyond the scope of this application note. 
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Figure 4: EMP FOURIER circuit 
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